;******************************************************************
FUNCTION mjd2yr, timecds,  decyr

;	1-Sep-2006	/gm


mjd=timecds.mjd

; 	find the year first, month, day first
	MJD2DATE, MJD, YEAR, MONTH, DAY
	
return, year
end

;******************************************************************
FUNCTION mjd2doy, timecds,  decyr


;	1-Sep-2006	/gm



	mjd=timecds.mjd      ; modified julian day

; 	find the year first, month, day first
	MJD2DATE, MJD, YEAR, MONTH, DAY
;	find the day of year
	
return, mjd-date2mjd(year,1,1)+1
end

;******************************************************************
FUNCTION mjd2decdoy, timecds,  decydoy 

;	1-Sep-2006	/gm


;	find the day of year
	doy = mjd2doy(timecds)

;	now build up fraction of the day


return, doy+timeCDS.time/(86400.*1000.)
end

;******************************************************************
FUNCTION mjd2decyr, timecds,  decyr


;	1-Sep-2006	/gm




mjd=timecds.mjd      ; modified julian day

; 	find the year first, month, day first
	MJD2DATE, MJD, YEAR, MONTH, DAY
;	find the day of year
	doy = mjd-date2mjd(year,1,1)+1
	if((year mod 4) eq 0) then ndays = 366.   else ndays = 365.

	
return, year+(doy + timecds.time/(86400*1000.))/ndays
end